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Abstract. We carry out an exact quantization of a PT symmetric (reversible) 
Lienard type one dimensional nonlinear oscillator both semiclassically and quantum 
mechanically. The associated time independent classical Hamiltonian is of non- 
' standard type and is invariant under a combined coordinate reflection and time reversal 

transformation. We use von Roos symmetric ordering procedure to write down the 
appropriate quantum Hamiltonian. While the quantum problem cannot be tackled in 
coordinate space, we show how the problem can be successfully solved in momentum 
space by solving the underlying Schrodinger equation therein. We obtain explicitly the 
eigenvalues and eigenfunctions (in momentum space) and deduce the remarkable result 
that the spectrum agrees exactly with that of the linear harmonic oscillator, which is 
also confirmed by a semiclassical modified Bohr-Sommerfeld quantization rule, while 
the eigenfunctions are completely different. 
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1. Introduction 

In a previous paper [Ij, Chandrasekar and two of the present authors have presented 
a conservative description for the Lienard type one dimensional nonlinear oscillator, 
namely 

X + kxx H x^ + uj'^x = 0, (1) 

9 

where overdot denotes differentiation with respect to t and k and u"^ are real parameters. 
Expressing ([T]) as a sysytem of first order equations, x = y = Fi{x,y), y = 
—kxy — ^x^ — uj'^x = F2{x,y), one can note that the divergence of the flow function 
F = fJ+ F2J of is non-zero (V.F = ^ + ^ = -kx). Also equation 1^ is 
invariant under the PT or reversible transformation, x — —x and t — >■ — t |2j. In spite 
of these, system admits a time independent Hamiltonian of the form 



H{x,p) 



2P 



2k f 2k Y k^x^ f 2A; 

-p - 2 1 - —p + — — 1 - —p 



3a;2 V 3a;2 / Ow^ \ 3a; 



-°°<P<^' (2) 



which is of non-standard type, that is the coordinates and potentials are mixed so that 
the Hamiltonian cannot be written as just the sum of the kinetic and potential energy 
terms alone, including velocity dependent terms. The corresponding Lagrangian L is 
given by 

_ ( 1 ^ 3a;2 . _ 9a;^ . 

^ - ^ (^A;x + fx2 + 3a;2 J + " 2F ' 

and the conjugate momentum is 

dL _ 21uj^ 3a;2 

9i ~ ^ 2k{kx + '^x^ + 3uj^Y ^ ~2k 

We note here that the system ([T]) also admits an alternate Lagrangian/ Hamiltonian 
(see for example Ref. [1]). However we consider the Hamiltonian given in the form ([2]) 
only since as — t- the Hamiltonian ([2]) reduces to the linear harmonic oscillator 
Hamiltonian, as the equation ([T]) does. We mention here that the parameters k and a;^ 
can be rescaled with appropriate scaling in x and t. However, to describe the physical 
properties of this system in the classical, semi-classical and quantum levels, we retain 
the parameters k and u"^ and do not scale them away. We also note here that the 
Hamiltonian (|2]) with the definition of p given by (jl]) is also invariant under combined 
action of coordinate reflection and time reversal (PT), x — )■ —x and t —t. 

The nonlinear oscillator (|2]) admits general periodic solution of the form 

, , Asm(ut + 6) „ , 3a; 

x(t) = rx^— 0< A<—, (5) 

^ ' 1 - If cos(a;t + 6) " ' ^ ^ 

where A, 6 are arbitrary constants. Note that for < A < ^ (while —00 < x < 00), the 
system ([1]) admits isochronous oscillations of frequency u, which is the same as that of 
the linear harmonic oscillator. For A> the solution becomes singular whenever the 
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phase (ut + S) takes the value cos~^ (|^) +2mi, n : any integer, even though x(t) given 
by the function in ([5]) is periodic of period T = while the corresponding momentum 
p is bounded and periodic (see equation ([7]) below). For more details on the classical 
dynamics of this system one may refer to Ref. [Ij- 

Exactly solvable quantum mechanical problems, particularly the ones involving 
nonlinear potentials are rare, even in one dimension. The few examples include 
Poschl- Teller, Morse, Scarf and isotonic oscillator potentials [3]. Also there exists a 
few velocity dependent potentials, for example Mathews-Lakshmanan oscillator and 
its generalizations [11[5]. The quantization of ([1]) is a challenging problem since 
the obstacles in this task are many. For example, the quantization of the damped 
linear harmonic oscillator itself is a quite complicated procedure requiring a rigged 
Hilbert space |6] description whereas the system under consideration is a nonlinear 
one. In addition to this, the associated time independent Hamiltonian is a non- 
standard one [7]. To the authors' knowledge there exists no nonstandard Hamiltonian 
system which is quantum mechanically exactly solvable. Such a system cannot also be 
quantized using standard techniques of canonical quantization P4T0]. In this paper, 
we completely solve the quantum mechanical problem of the Lienard type oscillator 
(II]), possessing the nonstandard Hamiltonian structure ([2]), by associating it with a 
position dependent mass Hamiltonian where now the variables x and p are interchanged. 
We then consider a general symmetric ordered form of the Hamiltonian proposed by 
von Roos [H] and solve the underlying Schrodinger equation in the momentum space. 
Since allowable choices of symmetric ordering lead to singular/unbounded solution, we 
transform the symmetric ordered Hamiltonian suitably in such a way that the associated 
Schrodinger equation possesses acceptable eigenfunctions. It is worth noting that the 
transformed Hamiltonian is now a non-symmetric ordered one as well as non-Hermitian. 
But it admits a real energy spectrum since the Hamiltonian and the corresponding 
eigenfunctions are invariant under the action of PT operation when — oo < p < Our 
results reveal that the eigenvalues of ([2]) exactly match with that of the linear harmonic 
oscillator, though the eigenfunctions are of a more complicated nature in the momentum 
space. The explicit form of the eigenfunctions is also presented. We also obtain the 
energy level spectrum through a semiclassical modified Bohr-Sommerfeld quantization 
rule for the regular periodic solution (|5]) which agrees with the quantum mechanical 
results. Additionally, we point out the existence of a negative energy spectrum in the 
quantum case corresponding to the sector p > 

We also note here the interesting fact that while the standard PT symmetric 
systems considered extensively in the recent literature [T2HT3] all correspond to PT 
invariant complex potentials involving complex valued dynamical variables, the present 
Hamiltonian system (|2]) is PT symmetric and real where the dynamical variables are also 
real. The motivation here is more of exact quantization of a nonlinear dynamical system. 
Consequently the analysis of the corresponding quantum system in the momentum space 
discussed below will also be different in spirit from the modified normalization scheme 
of complexified PT-symmetric schemes [15] . 
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2. Semiclassical quantization 



To start with let us consider ttie semiclassical aspects. To quantize the system 
semiclassically, we use the modified Bohr-Sommerfeld quantization rule [16], namely 



p dx 



(6) 



where h is the Plank's constant and n is any nonnegative integer and the integration is 
carried out over a closed orbit in the (a;,p) space. 

Firstly, we determine the energy of the system E = H using the general solution 
([5]) in ([2]). Plugging the expression ([5]) in (jl]), we obtain 



kA 

p = Au cos(ut + 6) { 1 cos(ut + 5) 



(7) 



Note that for regular (non-singular bounded) periodic oscillations in x{t) given by 
equation ([5]), the amplitude is restricted to the range < A < ^, while the range 
of p is restricted to — ^ < p < As we mentioned earlier (below ([5])) when A> 
one has singular periodic solution and there is no lower bound on p. Substituting now 
the expressions ([5]) and ([7]) in ([2]), we find that the energy of the system turns out to be 



E 



-AW. 
2 



Now, to evaluate the integral ([6]), we use ([7]) for p, and express dx from 



form 



dx 



Acos( 



out + S, 



so as to obtain the quantization condition for the regular periodic orbits. 



uoA' 



2-iT 



(cos 0-1^ cos^ 



0) (COS0-M) 



kA 



COS 



{n + -)h. 



(8) 
in the 

(9) 

(10) 



A^u 



n 



111 



(1 3. 

Evaluating the above integral, we arrive at 

2; ' - k 

Finally, from equations (j8]) and ( fTTl) . we obtain the allowed energy eigenvalues with 
an appropriate upper bound N onn corresponding to the regular periodic orbits of ([2]) 

as 

r 



Er, 



n + 



n 



0,1,2, ...AT 



(12) 



which agrees with that of the linear harmonic oscillator for this part of the spectrum. 
The semiclassical approach motivates us to prove that the energy of the nonlinear 
oscillator ([2]) can be quantized exactly as given in equation ( fT2|) . In the following, we 
proceed to solve the time indepedent Schrodinger equation associated with the system 
([2]) analytically, not in the coordinate space but in the momentum space. 
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Figure 1. The phase trajectories of the Hamiltonian system (IT5|) with w = fc = 1 for 
various values oi E = H . 



3. Quantum exact solvability 

Next we observe that the classical Hamiltonian H{x^ p) given in ([2]) is of the non-standard 
type, that is, 

1 



H{x,p) = ^f{p)x' + U{p), (13) 



where 



The (x — p) phase space structure is shown schematically in figure [H Note the 

hill ^ r, ^ 



deformed nature of the bounded periodic orbits around the origin and — %Tr < p < ^'^ 



3uj± 
2k 



The remaining trajectories have only an upper bound at p 

Note that in the limit A; 0, f{p) uj'^,U{p) ^, so that 

" = T + ^ (15) 

as it should be. Now the first term in the Lienard oscillator Hamiltonian fll3p contains 
both the position and momentum variables while the second term turns out to be a 
function of momentum alone. To quantize the Hamiltonian of this nature, one has to 
adopt a suitable ordering procedure. Since x and p are non-commuting variables in 
the quantum case, one may consider different ways of ordering between x and f{p) in 
order to quantize this Hamiltonian. After performing a detailed analysis we find that 
the nonstandard classical Hamiltonian given in ([2]) can also be equivalently considered 
in the form 
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where 

r-iP) = ^2 (1 1 Uip) = I, ^1 - |ip - 1 1 . (17) 

Interestingly, this form is similar to a position dependent mass Hamiltonian, H = 

2 11—1—1 

2 m(x) +^(^)? discussed extensively recently [Tl ]IT714T9] but with an important difference 
that the variables x and p are now interchanged. Once this fact is recognized one can 
consider a general symmetric ordered form of the quantum Hamiltonian proposed by 
von Roos in order to quantize the position (but now actually momentum) dependent 
mass Schrodinger equation [llj as 

H{x,p) = - [m" {p)x'm^ {p)x'm'^ (p) + m'^ {p)xm^ {p)xm°' (p)^ + U(p), (18) 

where the parameters a, /3 and 7 which remain to be fixed have to satisfy the condition 
a + /3 + 7 = —1. Obviously now in f[T8|) x and p are linear Hermitian operators satisfying 
the commutation rule 

[x,p] = ih. (19) 

We note here that since the variables x and p are interchanged we solve the 

Schrodinger equation corresponding to the Hamiltonian (ITSl) in momentum space with 
d 

X = ih— and obtain the time independent one dimensional Schrodinger equation in the 
op 

form 

H{x,p)iIj{x,p) = EiIj{x,p), (20) 



or 
2m 



m 



l + /3\/m'2 m"\ a(a + /3 + l)m" 



^ + 2— + ^ — ^ 



+ U{p)ij = E2P, (21) 



where prime stands for differentiation with respect to p. Since we are looking for bound 
states (even when U{p) becomes complex for p > ^), we can choose 

^ = 0, for p>— (22) 

and concentrate on the region —00 < p < alone in this section. We will also impose 
the boundary conditions 

V^(-oo) = ^ = (23) 

on the eigenf unctions for continuity and boundedness. 

Substituting now the expression for m(p) from flT7|) and its derivatives in equation 
( I2TI) and simplifying the resultant expression, we arrive at 

2k 1 , AkMc^ + P + , - -Bp-^) , 

^ 3c.^ (1 -^pf^ 9.^ (1 - ^pY"^ h^u^k^ (1 - ^p) ^ 

= 0, -00 < p < — . (24) 



Equation fl24p can be further simplified by introducing a transformation 
—p 01 y = \ 1 — -—^p- Since — oo < p < we nave < w < oo. 



1 - 



The resulting simplification of ( 12^ yields 

d'^tfj , 1# , , 4a(a + /3 + l) 



y 



-a(y-l)^)^ = 0, 0<y<oo, 



(25) 



where we have defined . „ , „ -E = and , „ , , 

k^ 

In order to solve (1251) subject to the boundary conditions corresponding to ( l23l) . 
namely ?/'(— 00) = ip{0) = 0, we first note the admissible asymptotic behaviour, 
ipi^y) — )■ e~^(^ as ?/ — )■ 00. So we introduce another transformation, namely 

in fl25|) so that it can be rewritten as 



(26) 



d^ 



dy' 



2 - 2^/Ey + 



1 

yj dy 



+ 



0. (27) 



y y 

Equation f l27|) can now be transformed to the Hermite differential equation under the 
change of variables 

<P{y) = y-'^'x{z), z = a'/\y-l), (28) 

with the condition Aa^a + (3 + 1) = —j. The transformed equation turns out to be of 
the form 



d'^Y dy 
^ — 2z— + 

dz^ dz 



With the restriction of the constant 



E- ^ 



a 



x = o. 



(29) 



2n, n = 0, 1, 2, equation fj29|) becomes 
the standard differential equation for the Hermite polynomials, 

d'x 



2z^ + 2nx = 0, 
dz^ dz 



n 



0,1,2,... 



(30) 



where x = Hn{z) which are nothing but the Hermite polynomials [20]. Then the 
eigenfunctions and eigenvalues can be readily written down as 



1pn{p) = Nn- 



exp 



_9a/ 
' 2 fc2 



2k 



xH, 



3a;3/2 



2k 



y/hk 
< p < 00, 



2k 



-p- 1 



—00 < p < 



3c^ 
2k 



and 



1, 



(n + -)huj, 



n 



0,1,2, 



(31) 
(32) 
(33) 



8 



Here N„ are constants. 



Note that the above eigenfunction is singular at the boundary p = ^ due to the 
denominator term in the right hand side of the equation (1311) and so the eigenfunction 
becomes unbounded. To avoid this singularity in the eigenfunction, we modify the 
starting Hamiltonian, H{x,p), suitably and solve the associated Schrodinger equation. 
For this purpose we can rewrite the time independent Schrodinger equation Hip = Eijj 
as 

so that we have 

1 



H ^ = {m Hm ) $ = and ^ = m iIj, m{p) 



(34) 
(35) 



^, one can have bounded, continuous and single valued 



(1 - m ■ 

With the choice d < ^ 
wavefunction for the Hamiltonian H. 

As a simple choice we consider a specific set of the values of the ordering parameters 
a = 7 = — I and /3 = — | that satisfy the conditions a + /3 + 7 = —1 and 
4a(Q; + /3 + 1) = — i so that d = Hence the Hamiltonian H{x,p) in f lTSj) becomes 



H{x,p) = -[m ^^^{p)xm ^^'^{p)xm ^^^{p)\ +U{p), 



(36) 



which admits the solution as given in ( 13T|) . Now we transform the Hamiltonian, given 
in ( |36l) . to the form 

(37) 



H = —;=H^/rri = -\rn '^^^{p)xm ^^'^{p)xm^^^ 



m 



—UJ' 1 



2k 

3a;' 



dp^ 



1 




. (3J 



This Hamiltonian (138|) is invariant under PT symmetry [2] though it is nonsymmetric 
and non-Hermitian. The Schrodinger equation corresponding to the Hamiltonian, 
H{x,p) is 

2 



2/l I _ 



2u2 



K'k 



9a;^ 
2P 



2k 

3u' 



1 <^ = E<^, 



24c^^ (1 - Bp) 

Equation ( 139|1 is now solved again by following the above procedure. The resulting 
bound state solution turns out to be 

L 1/4 f 9uj3 



Nn{l 



Bp) ' \ ^2hW 



1 



2k 
3^2 



p 



0, 



xH„ 



\^ k 



2k 



p — 1 



^-Bp 



— oo < p < 



3uj± 
2k ' 



(40) 



P^ 2k^ 

and the corresponding energy eigenvalues continue to be 

1, 



Er, 



n 



0,1,2, 



(41) 



9 



One can observe that the solution (Hn|) is continuous, single valued and bounded in the 
entire region — oo < p < oo and satisfy the boundary conditions $(oo) = <l>(3A;/2a;^) = 
oo) = 0. Since the eigenf unctions $n(p) are bounded and continuous in the region 
— oo < p < ^ and are zero outside this region, they are also normalizable. The 
normalization constants Nn can be found using the integration 



2k 



Kip) K{p)dp = l. 



This can be evaluated as 



On evaluating 



Kip) K{p)dp + 

becomes 



2k 



Kip) Kip) dp = 1. 



)n— 1 



mi 



aio) 



where 



9{a) 



3ui^ 
2k 



Kip) Kip) dp. 



Hence the normalization constant is 



1/2 



We further note that in the limit /c — )■ 0, equation fHOj) reduces to 



Kip) 



1/2 



exp 



2huj 



P 



P 



(42) 



(43) 



(44) 



(45) 



(46) 



-oo < p < oo, (47) 



which matches with the bound state solution of the harmonic oscillator in accordance 
with its classical counter part (vide (fT5|) ). 

Finally one can also note that one can choose many number of possible non- 
symmetric Hamiltonian H in fl37j) for suitable choice of the set of parameters a, /3, 7 and 
d all of which lead to the same eigenvalue spectrum but different sets of eigenf unctions. 



4. The p > ^ sector: broken symmetry 

In addition to the above solutions, we can also identify a different set of solutions which 
is nonzero only in the regime p > ^ with a different set of boundary conditions than 
( l23l) . To realize this, we consider solutions with $ = 0, —00 < p < and look for 
acceptable solutions in the region ^ < p < 00, either bound states with the boundary 
condition $(^) = = $(00) or continuum states or both. We also note that in the 
case of the classical nonlinear oscillator ([T]) with the Hamiltonian ([2]), there exists no 
real solution for p > ^ due to the form of the conjugate momentum (jl]). 
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For the above purpose, we consider the Schrodinger equation (15^ corresponding to 



the Hamiltonian (l38l) . under the transformation y = y — 1, as 
ci2$ ld$ / ~ 3 



- -— + [ -E + — + a{iy -ly ] ^ = 0, < y < oo, (4^ 



dy'^ y dy \ Ay"^ 

ISw^ ~ 3^*0;^ 
where again we have defined -E = E^ and = a. Under the transformation 

= v^e~~^ "'^'^x(z) with 2; = a"'^/'^(y + z), we get 

rfM-) 2.^-i±^xW=0, (49) 



so that the bounded solution (as \z\ 00 or y —> 00) can be now expressed (compared 
to ([29])) as 

x{z) = Hn{z), E = ^E = -{2n + l)V^, n = 0,1,2,.... (50) 



Thus the second set of solution to the Schrodinger equation ( 1391) can be written as 



(Bp - i)'''-P (-^& {Bp-^ + ^ ^^Bp-\ 



2fc ' 

0, -oo<p<|f, 
with energy eigenvalues without a lower bound 

En = -{n+^)hu, n = 0,1,2,3,.... (52) 

In (|5T]) Mn is the normalization constant. Note that the eigenfunctions (|5T1) are no 
longer PT symmetric, even though the Hamiltonian fl38|) is PT symmetric, leading to 
a negative energy spectrum that is unbounded below. Such a broken symmetry is 
obviously a consequence of imposition of a different set of boundary conditions for the 
sector p > ^ than ( 12 3 p which is reminiscent of the situation in the case of the linear 
harmonic oscillator [T2|[2T] and the H = — x'^ oscillator [12] . 



5. Conclusion 

We have shown that the non-Hermitian Hamiltonian H{x,p) given by ( |37l) admits 
the bound state solutions, ^n{p) and real energy eigen spectrum. En- The energy 
eigenvalues En also match with the energy values obtained through a semiclassical 
approach corresponding to regular periodic orbits. It is interesting to observe that 
the quantum system (!36l) possesses the energy eigenvalues En which are same as that 
of the linear harmonic oscillator, though the eigenfunctions are quite different from 
that of the linear harmonic oscillator. Our analysis shows that the underlying Lienard 
type PT-invariant reversible nonlinear oscillator is exactly quantizable and leads to 
interesting class of eigenfunctions and energy spectrum. It is also possible to generalize 
the above results to more general class of Lienard type nonlinear oscillators [22] and 
coupled nonlinear oscillators [23], which will be taken up elsewhere. 
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